Measurement of K*(892)° and K'^ mesons in Al+Al collisions at 1.9^4 GeV 

X. Lopez/'B N. Herrniann,2 K. Piasecki,^'^ A. Andronic,^ V. Barret/ Z. Basrak,^ N. Bastid,^ 
M.L. Benabderrahmane,^ P. Buehler,^ M. Cargnelli,^ R. Caplar,^ P. Crochet,^ P. Dupieux,^ 
M. Dzelalija/ L. Fabbietti,^ I. Fijal-Kirejczyk,^ Z. Fodor,^ P. Gasik,^ I. Gasparic,^ Y. Grishkin/° 
O.N. Hartmann,*^ K.D. Hildenbrand,'' B. Hong," T.I. Kang," J. Kecskemeti,^ M. Kirejczyk,^ 
Y.J. Kim," M. Kis,4'5 P. Koczon,-* M. Korolija,^ R. Kotte,!^ A. Lebedev/o Y. Leifels," V. Manko," 
J. Marton,^ T. Matulewicz,^ M. Merschmeyer,^ W. Neubert,^^ D. Pelte,^ M. Petrovici," F. Rami/^ 
W. Reisdorf,4 M.S. Ryu," P. Schmidt,^ A. Schiittauf,^ Z. Seres,^ B. Sikora,^ K.S. Sim," V. Simion,! 
K. Siwek-Wilczynska,'^ V. Smolyankin,^° K. Suzuki,^ Z. Tyminski,'^ P. Wagner, E. Widmann,^ 
K. Wisniewski,^ Z.G. Xiao,^^ I. Yushmanov," X.Y. Zhang,!^ A. Zhilin,i° and J. Zmeskal*^ 

(FOPI Collaboration) 

P. Kienle*^'^ and T. Yamazaki^'^ 

^ Clermont Universite, Universite Blaise Pascal, CNRS/IN2P3, Laboratoire 
de Physique Corpusculaire, BP 10448, F-63000 Clermont-Ferrand, France 
^ Physikalisches Institut der Universitat Heidelberg, Heidelberg, Cermany 
'"'Institute of Experimental Physics, University of Warsaw, Warsaw, Poland 
^ Cesellschaft fiir Schwerionenforschung, Darmstadt, Cermany 
"'Ruder Boskovic Institute, Zagreb, Croatia 
^ Stefan-Meyer-Institut fiir Subatomare Physik, Osterreichische Akademie 
der Wissenschaften, Boltzmanngasse 3, A-1090, Wien, Austria 
''University of Split, Split, Croatia 
^Physik Department, Technische Universitat Miinchen, D-85748 Carching, Cermany 
^KFKI Research Institute for Particle and Nuclear Physics, Budapest, Hungary 
''"institute for Theoretical and Experimental Physics, Moscow, Russia 
"Korea University, Seoul, Korea 
''^Institut fiir Strahlenphysik, Forschungszentrum Dresden-Rossendorf, Dresden, Cermany 
'"' Kurchatov Institute, Moscow, Russia 
'^'Institute for Nuclear Physics and Engineering, Bucharest, Romania 
'^Institut Pluridisciplmaire Hubert Curien, In2P3/CNRS, and Universite de Strasbourg, Strasbourg, France 
'^'Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou, China 
'^Heavy-Ion Nuclear Physics Laboratory, RIKEN, Wako, Saitama 351-0198, Japan 

(Dated; June 11, 2010) 

New measurement of sub-threshold A'* (892)" and K^ production is presented. The experi- 
mental data complete the measurement of strange particles produced in Al-I-Al collisions at 1.9 A 
GeV measured with the FOPI detector at SIS/GSI. The K* {9,<d2f / K'^ yield ratio is found to be 
0.0315 ± 0.006(stat.) ± 0.012(syst.) and is in good agreement with the UrQMD model prediction. 
These measurements provide information on in-medium cross section of K^--k~ fusion which is the 
dominant process on sub-threshold if* (892)" production. 



PACS numbers: 25.75.-q, 25.75.Dw 

The FOPI Collaboration has performed a high statis- 
tics experiment to study strangeness production in 
Al-l-Al collisions at a beam kinetic energy of 1.9A GeV. 
The collected data sample is large enough to enable for 
the first time the reconstruction of deep sub-threshold 
resonance production of i^*(892)". This measurement 
extends the data to four strange particle species re- 
ported in [l| with two strange resonances the S(1385) 
and the if* (892)°, and two neutral strange particles, the 
A and the . Strange resonances such as if* (892)° and 
I](1385) are particularly interesting to probe the earlier 
stage of the collision and the subsequent evolution of the 
medium. Indeed, these resonances are produced at sub- 
threshold energy and their productions stemmed mainly 



from the fusion of pions and K+ (A) J^, the latest being 
produced at the beginning of the collision [3, 4]. In addi- 
tion, if* (892)" and E(1385) have short life time (4 and 
5 fm/c, respectively) [5[. The measurement of if* (892)° 
and A'° is detailed in the present letter and their yield 
ratio is compared to the predictions of the UrQMD trans- 
port model. 

The experiment was performed at the Heavy-Ion Syn- 
chrotron SIS of the GSI in Darmstadt by using an Al 
beam of kinetic energy of 1.9yl GeV on an Al target. 
The target thickness was 567 mg/cm^ and the average 
beam intensity was 8 • 10^ ions/s. The FOPI detector 
is an azimuthally symmetric apparatus made of several 
sub-detectors which provide charge and/or mass determi- 
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nation over nearly the full solid angle. The central part of 
the detector is placed in a super-conducting solenoid and 
consists of a Central Drift Chamber (CDC) surrounded 
by a plastic scintillators (Barrel) for time of flight mea- 
surements. The forward part is composed of two walls of 
plastic scintillators (PLAWA and ZDD) and another drift 
chamber (Helitron) placed inside the super-conducting 
solenoid. More details on the FOPI apparatus and its 
different sub-detectors can be found in [6|-|8[ . The target 
was placed at 40 cm upstream of its nominal position in 
order to cover the whole backward hemisphere |9| . Events 
are selected according to their centrality which is deter- 
mined with the charged particle multiplicity measured in 
the CDC and in the PLAWA. The results presented here 
correspond to central collisions {(Jgeo = 315 mb) repre- 
senting about 20% of the total geometrical cross section. 

For the present analysis, the if* (892)° resonance 
(called K*'^ in the following) was reconstructed from its 
decay products in the channel K~^tt~ (branching ratio 
66% [Hi). Pions are measured in the CDC (23° < Oiab < 
114°) and are identified by their mass determined by 
the correlation between magnetic rigidity and specific en- 
ergy loss. are identified by matching the CDC and 
the surrounding ToF-barrel. This restricts the measured 
phase space of to 32° < 9iab < 57°. More details on 
identification with the FOPI detector are provided 
in 13, III- Long-lived neutral strange particle like the 



(30 for K+ and 40 for tt" 
track reconstructions; 



) for a higher quality of 



Kg are identified via their charged decay particles in the 
CDC. A complete description of the Kg reconstruction 
method can be found in [9] . 

The measurement of K*'^ corresponds to deep sub- 
threshold production (800 MeV) since the threshold en- 
ergy to create K*'-' in elementary reaction is 2.75 GeV. 
Because of its strong decay into K^ and 7r~, the short 
life time of this resonance does not allow us to disen- 
tangle its decay vertex from the collision vertex. There- 
fore, the invariant mass analysis consists of correlating 
primary kaons and pions. The invariant mass plot of re- 
constructed K*'^ is shown in Fig. [TJ In order to enhance 
the signal-to-background ratio, a set of conditions was 
imposed: 
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(i) a cut on K^ momentum (p < 0.68 GeV/c 
applied to reduce the contamination from fast pions 
and protons; 

(ii) a cut on the tt~ transverse momentum (pj"^ > 
0.185 GeV/c) is used to reject pions which are spi- 
raling in the CDC; 

(iii) a cut on the distance of closest approach (DC A) be- 
tween a track and the collision vertex in the trans- 
verse plane (DCA;^ + < 1 cm) was applied to 
suppress decay products outside the target; 

(iv) a condition was imposed on the minimum number 
of activated wires forming a track (hits) in the CDC 



(v) the following cut was applied for the systematic 
study of the signal (see text for more details): 
160° < e*j^+ + 91. < 200°, e* being the angle be- 
tween the momentum vector of a decay product in 
the reference frame of K*^ and the momentum vec- 
tor of K*^ in the reference frame of the collision. 
Its relevance is described in [l2|. 

The combinatorial background is obtained with the 
event-mixing method. The two decay particles {K^ and 
7r~) are taken from two different events characterized by 
the same particle multiplicity in the CDC. In addition, 
the two events are aligned to the reaction plane in or- 
der to have the same reference system for both particles. 
The reaction plane is estimated event by event utilizing 
the standard transverse momentum procedure proposed 
in |13| . The combinatorial background was normalized 
on the left and right side of the nominal mass of the 
K*^ . For either sides no difference was found on the 
signal characteristics (mean, width and counts). Figure 
1 shows the normalization on the left side (hatched re- 
gion in upper panel of Fig. [T]). The shape of the result- 
ing mixed-event background describes the combinatorial 
background and is indicated by dashed lines in Fig. [T] 
(upper panel) where the statistical error bars are smaller 
than the line thicknesses. After background subtraction 
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FIG. 1: (Color online) Invariant mass spectra of K'^iv~ 
pairs. The solid histogram and the dashed lines denote 
the data and the scaled mixed-event background, respec- 
tively (upper panel). The lower panel shows the signal af- 
ter background subtraction. The following characteristics of 
the signal are shown: number of counts in the signal (S), 
signal-to-background ratio (S/B) and significance (SIGNIF = 
S/\/S + B). The parameters extracted from the fit to the 
data (mean mass value (MEAN) and the width (F)) are also 
reported. 

(lower panel of Fig. [T|), within the extracted width (47 
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MeV), about 5400 are reconstructed from the anal- 
ysis of 290 million events. The fitted function used to ex- 
tract the number of counts of the signal is a Relativistic 
Breit-Wigner function modified by a Boltzmann factor to 
account for modification of the populated phase space in 
an hadronic heat bath 
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The mean mass and width 
extracted from the fit are in a good agreement, within 
statistical errors, with the values reported by the Particle 
Data Group [1| . The losses due to decay, acceptance and 
reconstruction efhciency have to be corrected in order 
to extract the particle yield. These corrections are deter- 
mined by means of GEANT simulations modeling the full 
detector response. The phase space distribution of 
was generated using the Siemens- Rasmussen formula (l5l | 
which describes a radially expanding equilibrated source 
characterized by a temperature T and a radial expansion 
velocity (3. The choice of the values of these parameters 
(T — 90 MeV, (3 = and 0.3) is imposed by previous 
analysis measurements [H, [3] and the values of the radial 
flow are detailed below. The K*'^ particles are embed- 
ded into an heavy-ion event calculated with the IQMD 
model 17|. The resulting output is subject to the same 
reconstruction procedure as for the experimental data. 
The distributions of simulated and measured event for all 
relevant quantities were carefully compared and found to 
be in good agreement Finally, the reconstruction 

efficiency is determined by computing the ratio of recon- 
structed particles in the simulation to those initially em- 
bedded into the background events. The systematic error 
on the K*'^ yield was evaluated in two steps. First the 
reconstruction efficiency was estimated from simulated 
K*^ with two values of the radial flow velocity (/3 = 
and (3 = 0.3). Then the K*^ signals were reconstructed 
under different sets of conditions on the relevant quanti- 
ties used to perform the reconstruction, as for example, 
different cuts on and n" momentum. The use of 
a new variable in FOPI analysis like the 9* [12], varia- 
tions of histogram characteristics (bin-size) and different 
normalization areas used to fit the background were also 
taken into account for the systematic calculation. 

Kg were reconstructed in the tt+tt^ decay channel 
(branching ratio 69.2% [1]) from its four momenta at the 
intersection point of the two pious tracks (Fig. [2|). The 
combinatorial background was evaluated, as for the K**^ 
case, with the event mixing method. After background 
subtraction, a signal of about 10^ counts within a ±2cr 
window is extracted. The large statistics of reconstructed 
candidates allows to perform a bi-dimensional {pt, y") 
efficiency correction in order to extract the yield of K'^, 
where is the reduced rapidity (?/° = (?/'"'' — 2/™)/y'^™). 
The procedure used to simulate Kg signals and Al+Al 
collisions is the same as the one used for the if*" par- 
ticle. Decay products of if ° cannot be measured with 
the FOPI apparatus. Since Kg and if ° are produced in 
equal amount (neglecting CP violation), the yield of if° 
corresponds to two times the one of Kg. The measured 
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FIG. 2: (Color online) Invariant mass distribution of 7r"'"7r~ 
pairs. The solid histogram and the dashed lines denote the 
data and the scaled mixed-event background, respectively 
(upper panel). The lower panel shows the signal after back- 
ground subtraction. See caption of Fig. [1] for the abbrevia- 
tions. 



transverse mass spectra (mt = ^/pf+rn^) of Kg after 
efficiency correction are shown on the left side of Fig . [S] 
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FIG. 3: (Color online) Left: Transverse mass distributions of 
Kg mesons for different rapidity windows after efficiency cor- 
rections. The lines represent the Boltzmann function fits. The 
range of the data used to perform the adjustment are shown 
by full symbols. The values in parentheses denote scaling fac- 
tors used for the plot. Right: Rapidity density distribution 
for where the open symbols are data points reflected with 
respect to mid-rapidity. The curve represents a Gaussian fit 
from which the total yield is extracted (see text for more de- 
tails). 

In order to extract the yield, we fit these distributions 
with a Boltzmann like function: 



1 



ml dirrit — mo)dy 



- = A ■ exp- 



-{mt - Too) 
Tb 



within a narrow window of rapidity dy'^. One can no- 
tice that, in the logarithmic representation, the measured 
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transverse mass distributions exhibit a linear decrease 
with increasing rrit — TOq and are described reasonably 
well by the Boltzmann function. The latter allows to 
extract the rapidity density distributions dN/dy^ by in- 
tegrating the fitted function from rrit — mo = to oo. 
The symmetry of the colliding system allows us to reflect 
the rapidity spectrum around mid-rapidity. The result, 
displayed in the right panel of Fig. [3l is fitted with a 
Gaussian function after applying correction factors by 
taking into account the branching ratio of (x 100/69) 
and the unmeasured (x2). The integration of the 
Gaussian distribution allows to extract the yield per 
central (315 mb) collisions: V(K°) = 0.038 ± 0.001(stat.) 
± 0.004 (syst.). The systematic errors of the measure- 
ment are obtained by using the same procedure applied 
for the K*° particles. Finally, the K*° to K° ratio is 
found to be: 



F{K 



*0\ 



P(/^o) 



0.0315 ± 0.006(stat.) ± 0.012(syst. 



Since for the K*^ resonance the efficiency correction is 
applied with a global factor due to the low statistics of 
the signal, such a ratio presents the advantage to cancel 
out in some extent the influence of the momentum dis- 
tribution chosen for the simulated signals. The depen- 
dence on rapidity of the Boltzmann temperature (Tg) of 
Kg is found to be in good agreement with an isotropi- 
cally emitting source (Tg = T^f f / cosh{yc.m.)) 0, where 
Ycm. is the rapidity in the nucleus-nucleus center of mass. 
For the mid-rapidity region, the inverse slope reaches 
Teff = 89 ± l(stat.) MeV and is consistent with the 
choice of T = 90 MeV used to generate simulated sig- 
nals. 

These measurements complete the study of strangeness 
production in Al-I-Al collisions at 1.9 A GeV as presented 
in [H. A summary of our experimental results is shown 
in Tab.lll The measured yield ratios are compared to the 
predictions by the UrQMD transport model 
The UrQMD model allows to simulate heavy ion reac- 



Yield ratios 


Data 


UrQMD 




5.3 ± 2.4 


6.9 


(A-hE0)/7r- 


0.017 ± 0.007 


0.018 


7^V(A + E°) 


0.66 ± 0.13 


0.78 


(S*+ +E-)/(A-HE0) 


0.125 ± 0.042 


0.177 




0.032 ± 0.013 


0.024 



TABLE I: Experimental ratios and predictions from UrQMD 
model. The error of the experimental results corresponds to 
the quadratic sum of statistical and systematic errors. 

tions according to a N-body theory on an event by event 
basis. It is the only QMD-based model which allows the 
production of S* and K* at SIS energy HJ. In this 



model, the Fermi motion (which is an important ingredi- 
ent for particle production below threshold) is treated in 



the following way. The initial momenta of the nucleons 
are randomly chosen between and the local Thomas- 
Fermi momentum The parameters used for the lat- 
ter are set to their default values as reported in [2^ . The 
model predictions are found to be in good agreement with 
the experimental ratios and their uncertainties. It is im- 
portant to note that there is no explicit in-medium modi- 
fications of particle properties in this transport code. The 
good agreement deduced from the comparison of model 
predictions and data could be explained by a low nuclear 
density reached during the collision in the light Al-f Al 
system. A detailed discussion concerning the E* reso- 
nance production is provided in [l|. The dominant K*^ 
production process involved in the UrQMD model is the 
fusion of kaons and pions (70%). The other processes to 
create K*'^ involve reactions of non-strange resonances 
(A and N*) with baryons. The time dependence of these 
reactions exhibits a maximum around 7.5 fm/c. This co- 
incides with the production time of as predicted in 
0. This imply that the kaon-pion fusion occurs directly 
after kaon production. 

In summary, we have presented new results on the 
P{K*'^)/P{K'^) strange mesons ratio in Al+Al collisions 
at a beam energy of 1.9 A GeV. For the first time, the 
strange resonance K*^ was measured 800 MeV below its 
production threshold in elementary reactions. The trans- 
verse mass spectra of Kg were found to be consistent 
with a Boltzmann distribution. The experimental results 
on strange particles measured with the FOPI apparatus 
were compared to the predictions of the UrQMD trans- 
port model. The measurements are in good agreement 
with the model predictions where the dominant process 
to produce if*" at sub-threshold energy is the fusion of 
kaons and pions. The measurement of the K* strange 
resonance in an heavier systems at sub-threshold energy 
could bring new informations on in-medium modification 
of particle properties at higher baryon density that can 
be reached at SIS energies and at the future FAIR facility. 
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